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Abstract 

We study the effects of bulk viscosity on pj spectra and elliptic flow in heavy ion collisions at 
RHIC. We argue that direct effect of the bulk viscosity on the evolution of the velocity field is 
small, but corrections to the freezeout distributions can be significant. These effects are domi- 
nated by chemical non-equilibration in the hadronic phase. We show that a non-zero bulk vis- 
cosity in the range (/s < 0.05 improves the description of spectra and flow at RHIC. 



1. Introduction 

The observation of nearly perfect hydrodynamic flow is one of the central discoveries of 
the heavy program at RHIC, and a significant amount of effort is being devoted to a precise 
determination of the shear viscosity to entropy density ratio rj/s. In this contribution we will try 
to estimate the bulk viscosity ( of the excited matter created at RHIC. 

Bulk viscosity enters the equations of fluid dynamics as an additional contribution to the 
stress tensor, ST*™ = -A ,iv (di l u k . Here, A /iv = gf v + u fl u v is a projector on the fluid rest frame 
and is the velocity of the fluid. Comparing with the stress tensor of an ideal fluid, T 1 ™ = 
(e + P)u M u v + Pg^ v where e is the energy density and P is the pressure, we observe that bulk 
viscosity reduces the pressure of an expanding fluid relative to its equilibrium value. In a heavy 
ion collision this implies that bulk viscosity reduces the amount of radial flow. 

Bulk viscosity also effects the spectra of produced particles. Particle spectra are computed 
by matching the stress tensor across the freeze-out surface. This yields the standard Cooper-Frye 
formula 

dN 

'd^p ~ (2n 

where dN/cPp is the spectrum of produced particles, E p is the single particle energy, and cr is the 
freeze-out surface. The distribution function f(E p ) = fo(E p ) + 6f(E p ) contains an equilibrium 
part /o and a viscous correction 6f. In the case of bulk viscosity 6f is proportional to the expan- 
sion rate duu k , but the overall magnitude and dependence on energy is sensitive to the underlying 
non-equilibrium reactions. 



E P • T^T f f&pWda-n , ( 1 ) 



2. Theories and models of the single particle spectra 

The bulk viscosity only constrains a moment of Sf, and determining the full functional form 
of the non-equilibrium distribution function requires a microscopic model or theory. The simplest 
model is based on the Boltzmann equation in the relaxation time approximation. In this approx- 
imation the complicated collision term in the Boltzmann equation is parameterized in terms of a 
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Figure 1: Non-equilibrium distribution^ of quarks and gluons in leading order perturbative QCD. The quantity^ is 
defined by Sf = -/o(l ± fo)x(p)(d ' «)■ The dashed curves show an approximate solution that does not exactly conserve 
energy. 



single collision time t(E p ). Energy and momentum conservation restrict the functional form of 
jjljE/;) llll- The bulk viscosity of an ultra-relativistic gas in the relaxation time approximation is 



( =l5\l-c 2 s ) V , (2) 



\2 

where c s is the speed of sound. The non-equilibrium distribution function is of the form 



We observe that the bulk viscosity scales as the second power of the conformal breaking param- 
eter (c 2 - 1/3), whereas Sf scales as the first power. This implies that for a nearly conformal 
fluid corrections to the freeze-out distribution are typically more important than corrections to 
the velocity fields. This conclusion does not depend on taking the relativistic limit. In general, 
the conformal breaking parameter is 

r = /^(T- c ^) /oW(1±/oW) (4) 

where + corresponds to bosons/fermions. The factor T vanishes in both the relativistic limit 
E p ~ p, and in the non-relativistic limit E p ~ m + p 2 /(2m). We find that £ ~ T 1 and Sf ~ T . 

The off-equilibrium distribution can be studied more rigorously in perturbative QCD j^. In 
QCD the process of emitting an extra soft gluon is efficient, and bulk viscosity is determined by 
the time scale for equilibrating the momenta of produced gluons via elastic 2 «-» 2 scattering. 
The off-equilibrium distribution of gluons with momenta p » T is 

^^r(H) (<9 - M) ' (5) 

where p^ ~ g 2 m 2 :) log(T /mo) is the drag coefficient and mo is the Debye screening mass. In 
perturbative QCD (c 2 s - 1/3) = 0(a 2 ). In pure gauge theory £ = 0A4a 2 T 3 / logia^ 1 ) which 




Figure 2: Transverse momentum spectrum of pions and protons at RHIC (data from | 6]). The lines show hydrodynamic 
calculations with and without bulk viscosity. 



implies f =2 48(c^ - 1 /3) 2 77. In full QCD we find interesting differences between the distribution 
functions of quarks and gluons, see Fig. Q] Both off-equilibrium distribution functions change 
sign, as is required by energy conservation, but the zero crossing occurs for different momenta. 

The difference between quarks and gluon distribution functions may manifest itself in vis- 
cous corrections to penetrating probes, but more direct observables are related to the differences 
between hadronic distribution functions. A problem that can be studied rigorously is the bulk 
viscosity of a pion gas |4J]. In this case bulk viscosity is determined by particle number changing 
processes, in particular the rate for the process n + n «-» An. In an expanding gas of massive pions 
the equilibrium value of the total pion number decreases with time, but pion number changing 
processes are slow and a pion excess is created. The distribution function can be parametrized 
by an off-equilibrium chemical potential for the total number of pions 

sf = fo ( y + ^jP] = -f° (*° ~^ e p) ^- u ^ ^ 

where 6fi is related to the bulk viscosity f and ST is fixed by energy conservation. The bulk 
viscosity is controlled by the inelasticity AE = 2m K , £ ~ (f% /ml)exp(-2m„/T). We have ex- 
tended equ. (O to a hadronic resonance gas, see for earlier studies of chemical equilibration. 
We assume that the relative sizes of dp for different species are determined by fast reactions like 
p <-> nn or p + p <-> nn (n ^ 5 in the regime of interest). This implies, for example, fi p = 2p n and 
p p = 2.5fi n . The overall scale of 6/x is related to the bulk viscosity £ and can be extracted from 
experiment. As before 6T is determined by energy conservation. 

3. Spectra and flow at RHIC 

We have applied this model to pj spectra and flow at RHIC, see Figs.|2]and[3] The hydro- 
dynamic model incorporates a single freezeout temperature 7y> 140 MeV and no hadronic 
afterburner, but the spectra include feed-down from hadronic resonances. The model for the 

3 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 3: Differential elliptic flow of pions and protons at RHIC (data from [7]) compared to hydrodynamic calculations 
with and without bulk viscosity. 

temperature dependence of the bulk viscosity is described in U. The value of £/ s quoted in the 
Figure refers to the freezeout surface. We have solved second order hydrodynamic equations and 
we have verified that the gradient expansion is convergent. Viscous corrections to the spectra 
become large for pj > 2 GeV, and results in this regime cannot be trusted. We observe that 
a non-zero bulk viscosity improves the description of the spectra and flow. In particular, bulk 
viscosity raises the single particle spectra at low pj, and increases the splitting between the pion 
and proton Vzipr)- These effects cannot be described in purely hydrodynamic models without 
bulk viscosity, but they have been explained in terms of hadronic non-equilibrium effects in ki- 
netic afterburners. Our results show that these effects can be described efficiently in terms of 
bulk viscosity. We also note that many important hadronic reactions, such as pp annihilation 
into several pions, are difficult to include in kinetic models. Finally, we observe that the value 
of £ extracted from the data is surprisingly small, £/s - 0.01. The corresponding pion chemical 
potential is of the order p n (10 - 20) MeV, depending on the local expansion rate. 
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